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Abstract 

Linear and angular momenta of a soliton in a ferromagnet are commonly derived through 
the application of Noether’s theorem. We show that these quantities exhibit unphysical 
behavior: they depend on the choice of a gauge potential in the spin Lagrangian and can 
be made arbitrary. To resolve this problem, we exploit a similarity between the dynamics 
of a ferromagnetic soliton and that of a charged particle in a magnetic field. For the 
latter, canonical momentum is also gauge-dependent and thus unphysical; the physical 
momentum is the generator of magnetic translations, a symmetry combining physical 
translations with gauge transformations. We use this analogy to unambiguously define 
conserved momenta for ferromagnetic solitons. General considerations are illustrated on 
simple models of a domain wall in a ferromagnetic chain and of a vortex in a thin film. 


1. Introduction 

The definition of linear and angular momenta of a ferromagnet is a subject with a 
long history, surprising results, and a lingering controversy [1—12]. It was realized in 
the 1970s that the linear momentum of a ferromagnetic soliton is determined not by its 
velocity but rather by its configuration [1, 2]. Although this sounds counterintuitive, one 
needs to realize that our physical intuition is based on experience with massive objects, 
for which an external force generates proportional acceleration. Spins in a ferromagnet 
behave differently: like fast-spinning gyroscopes, they precess at an angular velocity 
proportional to the external torque. An external force F acting during a short time 
interval dt increments the velocity of a Newtonian particle, Fdt = mdv. For a spin, 
an external torque t affects its orientation, rdt = dS. Thus the linear momentum of a 
ferromagnetic soliton is a function of its collective coordinates, rather than velocities. 

The earliest derivations of linear momentum in a ferromagnet followed the above qual¬ 
itative reasoning and analyzed the configurational change of a soliton under a specified 
external perturbation. For example, a domain wall in a ferromagnetic chain pushed by 
an external magnetic field increments its azimuthal angle in proportion to the impulse 
of the force exerted by the field. The linear momentum of the domain wall is therefore 
proportional to its azimuthal angle [1]. 

On a deeper level, momenta are conserved quantities related to global symmetries. 
By Noether’s theorem, invariance of the Lagrangian under translations and rotations 
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Figure 1: Magnetic field b of a magnetic monopole with a net flux —4irJ (blue arrows). A Dirac string 
at m s (tube) carries an equal and opposite magnetic flux 4-4x L 7 (red arrow). Standard choices are 

Ills = ±Z. 


gives rise to the conservation of linear and angular momenta. The precessional dynamics 
of the magnetization held m(r, t) of unit length is represented in the Lagrangian not by 
a kinetic energy, but rather by a Berry-phase term Lb = a ■ <9 t m. Here a(m) is the gauge 
potential of a magnetic monopole [13] whose held is 

b(m) = V m x a(m) = -J\ n. (1) 

J is spin density. Eq. (1) does not fully specify a(m): any gauge transformation 

a'(m) = a(m) + V m y(m) (2) 

preserves b = V m x a. It is known that the linear momentum p derived from Noether’s 
theorem is gauge-dependent [8, 11]. This is a cause for concern: physical quantities 
should be gauge-invariant. 

A different but related problem arises for angular momentum J. Although the mag¬ 
netic held of a monopole (1) is spherically symmetric, its gauge potential is not. The 
best we can do is to make it axially symmetric [3]: 

L b = a • d t m = J- - • d t m. (3) 

1 — m, • m 

The singular direction m s is the location of a Dirac string carrying away the magnetic 
hux 471^7, Fig. 1. The axial symmetry of the Berry-phase term in the Lagrangian limits 
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us to just one conserved component of angular momentum J ■ m s . The string is usually 
placed at m s = ±z, which yields the familiar kinetic term [13] 

C B = a • d t m = J(cose ± 1) d t <t>. (4) 

We shall show below that linear and angular momenta defined through Noether’s the¬ 
orem for a Lagrangian with the kinetic term (3) depend on the location of the string m s . 
Momenta defined in this way are essentially arbitrary, and thus unphysical, quantities. 

Our goal will be to define conserved momenta that are free from these problems. 
The resolution is tied to a powerful analogy between the dynamics of a magnetic soliton 
and of a charged particle in a magnetic field pointed out by Bar’yakhtar and Ivanov 
and by Papanicolaou [6]. For a particle in a magnetic field, the gauge potential A 
breaks translational invariance even if the field B=VxA itself is uniform. Canonical 
momentum p = dL/di = mr + eA is neither gauge invariant, nor conserved. The proper 
symmetries of the Lagrangian are combinations of translations and gauge transformations 
known as magnetic translations [14]. Their generator 

P = p — eA — er x B = mi — er x B (5) 

is both gauge-invariant and conserved. A related quantity, the guiding center of the 
cyclotron orbit r 0 , is defined by the equation r = (r — r 0 ) x eB/m [15]. 

Conserved momentum of a particle P = mi — er x B has a transparent physical 
meaning. Newton’s second law contains three terms, 

mi = er x B + F, (6) 

where the Lorentz force er x B has been set aside. Treating it as a kinematic effect, just 
like the inertia term mr, gives the right perspective. The impulse of the external force 
F then defines the change in momentum. 

AP = f J F dt = A(mi — er x B). (7) 

Hence Eq. (5). 

Below we will first show that the momenta of a ferromagnetic soliton derived through 
the naive application of Noether’s theorem are canonical momenta conjugate to the 
soliton’s zero modes. They are neither gauge invariant, nor conserved, generally speaking. 
We will then derive the physical conserved momenta that are analogs of the generator 
of magnetic translations. The method is applicable to any magnetic soliton and is not 
restricted to the symmetries of translation and rotation. 

The paper is organized as follows. Sec. 2 introduces the general formalism of La¬ 
grangian mechanics for solitons in a ferromagnet and defines gauge-dependent canonical 
momenta and their gauge-invariant conserved counterparts. The formalism is illustrated 
on the examples of a domain wall in one dimension (Sec. 3) and of a vortex in two 
dimensions (Sec. 4). Sec. 5 summarizes the main results. 

2. Mechanics of a soliton 

2.1. Lagrangian and the gauge field 

A soliton can be described in terms of some collective coordinates { 91 , 92 , ■ • •} such 
as its position, size, orientation, shape, and so on [16]. Time evolution of the mag¬ 
netization field m(r, t) can then be expressed as the evolution of these coordinates: 
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c^m = q t dm/dq.i. Here g, = dqi/dt ; summation is implied over doubly repeated indices. 
The Lagrangian L = fa ■ d t mdV — /7[m], expressed in terms of collective coordinates, 
becomes [17] 

L = Ai(q)qi - U(q), A* = J dV a(m) ■ dm/dq.,. (8) 

The Lagrangian is similar to that of a massless particle moving in an external magnetic 
field, whose strength F^ is given by the curl of the gauge potential A, : 

Fij = diAj - djAi. (9) 

U is potential energy and di = d/dqi. 

To prevent confusion, let us note that the Lagrangian with no mass terms (8) applies 
when we keep all of the system’s collective coordinates, not just its soft modes. If we 
integrate out some hard modes, soft modes may acquire inertia. This is the origin of 
Doring’s mass [18, 19] observed in domain walls [20] and skyrmions [21]. 

To make contact with the more familiar case of a massive particle, we will add an 
inertia term, 

L = mqiqi/2 + Aiqi - U, (10) 

and will set m = 0 when we deal with magnetic solitons. 

2.2. Canonical momenta and conserved momenta 

Continuous symmetries are associated with cyclic coordinates, which do not affect 
the potential energy and the magnetic field: 

d a u = 0, d a Fij = 0. ( 11 ) 

To distinguish these cyclic coordinates from the rest, we reserve for them indices from the 
beginning of the alphabet, a, 6 ,... Canonical momenta conjugate to cyclic coordinates 
are 

p a = dL/dq a = mq a + A a . (12) 

They are neither gauge invariant, nor conserved: p a = d a Ai cp ^ 0 for a generic gauge 
choice. 

To construct conserved momenta, we examine the transformation of the Lagrangian 
under an infinitesimal translation: 

Sq a e a , dAj e a d a Ai, SL q^SA^. (13a) 

The change in the Lagrangian stems from the lack of translational invariance of the gauge 
potential, A(q + e) ^ A(q). However, the translated potential A(q + e) describes the 
same magnetic field, F(q + e) = F(q), by virtue of translational symmetry (11). We may 
then follow up the translation (13a) with a gauge transformation 

5Ai = dix( q), SL = & = d\/dt, (13b) 

to return the gauge field, and with it the Lagrangian, to its original form. (Adding a 
full time derivative to the Lagrangian does not affect the equations of motion.) Con¬ 
served momentum is the generator of the combined symmetry of translation and gauge 
transformation, 


P a = dL/dq a + dx/de a 
4 


( 14 ) 


To find the right gauge transformation (13b), we solve the equation di\ + e a d a Ai = 0 
for x- To that end, we rewrite d a Ai = F al + diA a and obtain a tentative answer 


x(q) 


-z a A a { q) - e-a 


F a i{ q) dq t . 


(15) 


It remains to verify that the second term in Eq. (15) depends on the initial and finite 
positions but not of the integration path between them. That would be the case if the 
integral over any closed loop vanishes, j> F ai dqi = 0. In the differential form, d 3 F ai — 
diF a j = 0 (zero curl). That this condition is fulfilled can be seen by noting that the field 
is uniform in the q a direction, d a Fij = 0, and that it satisfies the Jacobi identity 

diFjk + djFki + dk^ij = 0 (16) 

as long as it is derived from a gauge potential, Fij = diAj — djAi. Upon combining 
Eqs. (14) and (15), we obtain the conserved momentum 


P a (q, q) = p a - A a 


Fai dqi = mq a - 


Fai dqi . 


(17) 


This quantity is both gauge invariant and conserved. Eq. (5) is a particular case of a 
particle in 3 dimensions in a uniform magnetic field, F a b = e a b c B c . 

An alternative way to derive the conserved momentum is to start with Newton’s 
second law, 

mq a = F^+F^\ (18) 

where the right-hand side includes the Lorentz force F^ = F a iiji and the conservative 
force Fa^ = —d a U. For a cyclic coordinate, the conservative force vanishes. Let us break 
the symmetry and apply an external force F^ 7 ^ 0. Its impulse induces a change in the 
momentum P a . With the aid of Eq. (18), we obtain 

A P a = J F^ c) dt = mAq a - J F ai dqi, 


in agreement with Eq. (17). Thiele [2] followed similar logic to define linear momentum 
of a domain wall. 

For a soliton in a ferromagnet, the mass to = 0, so conserved momenta are functions 
of coordinates: 

r<i2 

Pa{ q2) - Pa{ qi) = - / Fai{ q) dqt. (19) 

•' qi 

The role of the magnetic field F t j is played by the gyrotropic tensor for collective coor¬ 
dinates [17] 


Fi, 


J 

J 


( dm 
\dqi 


/ d(j) d cos 6 
\ dqi dq 3 



d(j> <9cos0\ 

dq-j dqi J 


dV. 


( 20 ) 


This expression can also be obtained directly from the definitions of the gauge potential 
( 8 ) and magnetic field (9), see Appendix A. 
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It is useful to translate the definition of momentum (19) from the language of collective 
coordinates back to the field-theoretic description. To compute the difference of momenta 
P a between two arbitrary states nii(r) and m 2 (r) of a magnetic soliton, we imagine 
evolving it in time from the initial state m(r,fi) = mi(r) to the final state m(r,f 2 ) = 
ni 2 (r). We express F a i dqi = F a i qidt with the aid of Eq. (20) to obtain 

P a [m 2 (r)] - P a [mi(r)] 

The momentum difference depends only on the initial and final configurations but not 
on the path between them, nor on how fast the evolution happens. 

Canonical momentum is 


p a = dL/dq a = A, 


d V a(m) 


dm 

dqi ' 


2.3. Some useful relations 

The standard Poisson bracket is defined as 

r, i = ( df dg df dg \ 

U,9} ^ \ dpi dq t dqi d Pi ) ' 


( 22 ) 


(23) 


Poisson brackets for coordinates and conserved momenta of a massive particle are thus 
{<?*, Qj} = 0, { Pa , qi} = S ai , {P a , Pb} = ~F ab . (24) 

For magnetic solitons, the mass term vanishes so that canonical momenta (12) lose 
their dependence on velocities and become functions of coordinates only, p a = A a (q). 
Because coordinates and canonical momenta are no longer independent, the standard 
definition of the Poisson bracket (23) breaks down and needs to be modified. 

In a ferromagnet, the z component of a spin HS cos 9 is the canonical momentum 
conjugate to its azimuthal angle 4>. Hence the Poisson bracket of a discrete set of spins 
{S r } [22], 

df dg dg df 


{/,<?} = 

r 

and the continuum version, 


d cos 9 r d(p r d cos 9 r d<j> T 


= -j 


Sf Sg Sg Sf 


6 cos 8 5 cj) <5 cos 9 S(j> 


dV. 


(25) 


(26) 


The Poisson brackets for collective coordinates and conserved momenta of a ferromagnet 
are (Appendix B) 

{quqj} = {P a ,qi} = S a i, {P a ,Pb} = ~F ab - (27) 


Eq. (27) indicates that at the quantum level conserved momenta of a soliton do not 
commute if the corresponding gyrotropic coefficient is nonzero. Watanabe and Murayama 
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Figure 2: A domain wall in a ferromagnetic wire. 


[23] pointed this out for a ferromagnetic skyrmion. We see here that non-commutativity 
of conserved momenta in a ferromagnet is a generic feature. Furthermore, coordinates 
are non-commutative, too! 

The following gauge-invariant relations will be helpful: 

dpb _ dpa 
dq a dq b 
dP b dP a 
dq a dq b 

They follow directly from the definitions of conserved (19) and canonical (22) momenta 
of a soliton. Eqs. (28) indicate that it is impossible to find a gauge in which two canonical 
momenta are equal to their conserved counterparts, p a = P a and p b = P b , unless the 
gyrotropic coefficient F ab vanishes. 


Fab, 

(28a) 

2 F ab . 

(28b) 


3. Domain wall in one dimension 


3.1. Lagrangian and soliton solutions 

A ferromagnetic chain with easy-axis anisotropy has a Lagrangian density 


C = a(m) • <9 t m — 


A(d z m) 2 A'(m x z) 2 


(29) 


The system has two uniform ground states m(z) = ±z. We shall consider here topological 
solitons satisfying the following boundary conditions: 

m(±oo) = ±z. (30) 

In equilibrium, a domain wall has the width A = \JA/I\ and configuration 


m x + im y 



m z = tanh 


2- Z 

A 


(31) 


The two free parameters—the location of the wall Z and the azimuthal angle <1>—are 
zero modes associated with the symmetries of translation and spin rotation about the z 
axis. 
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Pz = 




Pz = 2t:J Pz = -27r 

Pz 



Figure 3: A one-dimensional domain wall in equilibrium (31) is depicted on the unit sphere of magneti¬ 
zation m as an arc connecting the south and north poles along the 4> meridian. With the Dirac string 
(tube) away from the poles, m s 7^ =Lz, canonical momentum pz (34a) increases linearly with <f> as the 
arc moves east, from pz = —2tt J just east of the string to +27r J just west of the string. 



















3.2. Canonical momenta 

Canonical linear momentum is 

/ +oo p+z 

a{in) ■ d z mdz = — / o(m) • dm, (32) 

-OO J —Z 

with the integration path between the poles taken along the meridian $ = const. In 
deriving this result from Eq. (22) we relied on the fact that d/dZ = —d/dz for a rigidly 
translated domain wall m (z—Z). The expression for canonical momentum (32) coincides 
with that derived through the application of Noether’s theorem [ 8 , 11]. It changes under 
a gauge transformation ( 2 ): 

/ +z 

V m x(m)-dm=p z -x(+z)+x(-z)- (33) 

Z 


We see that canonical momentum pz of a domain wall m(r) is not well defined. 

We may attempt to remove the unphysical gauge dependence by evaluating the dif¬ 
ference of momenta between two states of a domain wall mj(z) and mjfz) [ 8 , 11]. The 
gauge dependence cancels out, provided that the two configurations are in the same 
topological sector (30). 

However, we are not out of the woods yet. Although momentum difference pz[m 2 (2)] — 
pz[m 2 {z)\ does not change after a gauge transformation, this quantity is still sensitive 
to the location of the Dirac string m s in Eq. (3) as we shall see next. As the string is an 
artifact of the gauge description, the dependence on its location makes canonical linear 
momentum unphysical even after the subtraction. 

For an axially symmetric gauge (3) with the Dirac string away from the poles, we 
obtain 

$ — d> s 

pz = —AJ arctancot—-— 

= 2J[<& — <j) s — (fin + l) 7 r], m s 7 ^ ±z, (34a) 


for $ between <j> s + 2mr and 4> s + 2(n + l) 7 r. Momentum increases linearly with $ from 
pz = —2-kJ with the domain wall just east of the Dirac string to pz = +2 nj just west 
of the string, Fig. 3. Crossing the string results in a jump A pz = —AirJ. Thus the 
presence of the string ensures the single-valuedness of canonical momentum. 

The AmJ discontinuity of canonical momentum is a bug in the classical field theory 
of a ferromagnet. It goes away when we recall that the continuum theory is derived from 
a system with discrete quantum spins [3]. Translations become a discrete symmetry. In 
a chain with a lattice spacing a and spins of length S, the spin density is J = hS/a. The 
generator of translations is the translation by one lattice spacing, T(a) = exp (— ipza/h ). 
Adding Anff to pz multiplies the translation operator by 


exp (—Anija/K) = exp (—47 riS) = 1 


(34b) 


if S is integer or half-integer. We see that the discontinuity does not affect the physical 
operator T(a). 

Although the appeal to the discrete and quantum nature of the ferromagnet saves the 
day, the solution comes at a high price of abandoning the classical field theory. There 
are other problems, too. 
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When the Dirac string is located at one of the poles, m s , = ±z, Eq. (34a) does not 
apply as the azimuthal coordinate of the string <j) g is undefined. Evaluation of Eq. (32) 
in the gauge (4) yields 

pz =0, m s = ±z, (34c) 


because d(j) = ( dcj)/dz)dz = 0 in equilibrium, see Eq. (31). Linear momentum vanishes 
in this gauge. 

Canonical angular momentum is equally problematic. To evaluate it, we introduce 
a local reference frame with mutually orthogonal unit vectors pointing up, south, and 
east: 


<9m 

e ! = nr e 2 = —, e 3 


1 9m 
sin 6 d(f> 


(35) 


With the aid of Eq. (22) we obtain 


= -J 



m s ■ e 2 sin 9 
1 - m s • e x 


dz. 


For a domain wall in equilibrium (31), — e 2 sin 9 = A dei/dz, which allows us to evaluate 
the integral, 

= -JX ln[l - m s ■ m(z)]|+“ . (36) 

For a generic location of the Dirac string away from the poles, canonical angular 
momentum is a gauge-dependent constant: 


= 2J A In cot , 


m s ^ ±z. 


(37a) 


With the string at one of the poles, the result diverges and we need to work with a chain 
of large but finite length L A. Working in a standard gauge (4) yields 

r+L/2 

= J / (cos0±l)dz. 

J-L/2 


The same expression can be obtained via Noether’s theorem [7, 10]. For a long chain, 
L » \Z\, L » A, 

= ±JL — 2JZ, m s = ±z. (37b) 


Like its linear counterpart, canonical angular momentum shows strong gauge depen¬ 
dence. In the standard gauges (4), depends on position Z. For other positions of the 
Dirac string, it does not. 

The dependence of canonical momenta (34) and (37) on the location of the Dirac 
string clearly makes them unphysical. It is not obvious a priori which of the answers, if 
any, is correct. Our best bet is to avoid them entirely and to use conserved momenta 
that are free from these artifacts. 
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3.3. Conserved momenta 


Gauge-invariant conserved momenta Pz and for a domain wall in equilibrium 
can be computed with little effort. Because no other (hard) modes are excited, we only 
need the gyrotropic coefficients (20) involving both soft modes, F$z = —F z <s> = 2 J [17]. 
Using the domain wall with Z = $ = 0 as a reference point, we obtain with the aid of 
Eq. (19) 


Pz(Z,&) = -Pz( 0 , 0 ) + 2^7$, (38) 

P*(Z,$) = P*(0,0)-2JZ. (39) 

The proportionality of angular momentum to coordinate Z is easy to understand. In 
a ferromagnetic chain, angular momentum comes from spin alone. Shifting the domain 
wall from 0 to Z elongates the m = — z domain and shortens the m = +z domain by 
Z , thereby reducing the z component of spin by 2 JZ. This result is not sensitive to the 
detailed structure of the domain wall in equilibrium (31) and remains valid as long as the 
soliton interpolates between the two ground states, m(±oo) = ±z, and the azimuthal 
angle is spatially uniform, (j>(z) = $. 

Unlike its canonical counterpart, conserved momentum Pz (38) is not single-valued: 
advancing the azimuthal angle d> by 27 t increments Pz by FkJ even though the domain 
wall returns to the original state. The multi-valuedness is connected to the absence of 
the Dirac string in the gauge-invariant treatment and to the ensuing violation of the 
Jacobi identity (16). 

Comparing the canonical momenta obtained in various gauges (34) and (37) with 
the physical answers (38) and (39), we observe that canonical momenta p a sometimes 
reproduce the correct answers P a and sometimes they do not. Furthermore, no gauge 
choice yields pz = Pz and = P$. The canonical approach gets right at most one 
or the other, but not both. That is not a coincidence: Eqs. (28) show that, as long as 
Fab 7 ^ 0 , there is no gauge in which two canonical momenta p a and pb can match their 
physical counterparts P a and Pb- 

One may wonder how Yan et al. [11] managed to get both linear and angular momenta 
right in the same gauge. The answer is their linear momentum should have been 0 as in 
Eq. (34c). As explained by Thiele [2], the formula pz = J J (cos 6 ± 1) d<j> in a standard 
gauge (4) is missing crucial boundary terms. Yan et al. [11] smuggled them in. 

To compute the linear momentum of a domain wall with arbitrary deformations, we 
use Eq. (21) and replace d/dZ with —d/dz (as Z represents a rigid displacement) to 
obtain 


Pz [m 2 (z)\ - P z [mi (z)] = J JJ x ^dtdz = JA. (40) 

The double integral represents the area A swept by the curve m (z,t) on the unit sphere 
as it evolves from mi(z) to m 2 (z). This result was correctly anticipated by Galkina and 
Ivanov [8] and by Yan et al. [11]. If both the initial and final states are equilibrium 
configurations (31) then the momentum difference is 2j7(4>2 — 4>i), in agreement with 
Eq. (38). Conserved linear momentum for a domain wall (40) was first obtained by Thiele 
[2], who derived it by integrating the impulse of the gyrotropic force. 
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X = +l X = -l 




Figure 4: Four possible equilibrium configurations of a vortex in a thin film. Arrows represent the in¬ 
plane components m x and m y , their colors the azimuthal magnetization angle <j>, and their brightness 
the out-of-plane component m z (black is — 1 , white is +1). 


4. Vortex in a thin film 

4-1. Soliton solutions 

In a thin film, shape anisotropy due to dipolar interactions forces magnetization m 
into the plane of the film, m z = 0. A vortex possesses a small core where m z 7 ^ 0. 
Away from the core, the in-plane magnetization rotates through +27r as the core is 
circumnavigated counterclockwise. In equilibrium, a vortex has an axially symmetric 
structure. In polar coordinates (r, a) with r = 0 at the core, m z is a function of the 
radius r, whereas the azimuthal magnetization angle <f> is a function of the azimuthal 
angle in configuration space a: 

m z ( 0) = p = ±1, 777 - 2 ( 00 ) = 0, (41a) 

4>(a) = a + x'7t/ 2 = a ± 7 t/ 2, (41b) 

where \ = ±1 is the “chirality” of the vortex and p = ±1 the “polarity” of its core. 

The locking of magnetization to the azimuthal direction in configuration space (41b), is 

enforced by dipolar interactions [24]. Two values of polarity and two of chirality give 
four distinct vortex configurations, Fig. 4. 
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J^.2. Canonical momenta 

Canonical momenta can be obtained from Eq. (22) along the lines of Sec. 3.2. Linear 
momentum is 

p a = —J / a(m) • d a md 2 r. (42) 

J n 

Here H is the two-dimensional area of the sample. Under a gauge transformation (2), 
linear momentum changes as follows: 

P' = P- / V r xd 2 r = p - / xndr. (43) 

Jn Jd n 

Here dfl is the boundary of the sample and n is its outward normal. Again, canonical 
momentum is gauge-dependent and unphysical. 

In the standard gauges (4), 


p a = —J / (cos 6 ± l)d a (j)d 2 r. 


(44) 


The same expression results from the application of Noether’s theorem [5, 6]. Evaluation 
of the integral in Eq. (44) requires care: the integrand decays slowly away from the core. 
To ensure its convergence, we restrict the integration to a disk of radius Rd much larger 
than both the vortex core and the vortex displacement R and subtract the momentum 
of the vortex at the origin, p a (R) — p a (0) (see Appendix C). Upon the subtraction, con¬ 
tributions from the interior of the disk cancel out, leaving behind an integral over a shell 
of width R <C Rd on the edge. The difference is linear in R and is independent of vortex 
polarity p because the contributions from the core canceled out after the subtraction: 

p a ( R) = Pa(0) ± 7 rJeapXp, m s = ±z. (45) 

To avoid the spurious contribution of the Dirac string to canonical momentum p a , 
we should attach the string where magnetization cannot reach it, i.e., opposite to the 
magnetization at the core, m s = —pi. Then 

Pa (R) = Pa (0) - pnje a pXp. (46) 


We shall see shortly that these values are only half the correct answer, given below by 
Eq. (52). As explained in Sec. 2.3, this failure is unavoidable. That canonical momentum 
p is only half of conserved momentum P is consistent with Eqs. (28). 

Like with the domain wall, it is possible to match one of the canonical momenta with 
its conserved counterpart—at the expense of the other. Place the Dirac string somewhere 
at the 4> s = 0 meridian, m s = (sin0 s ,O,cos0 s ). We then obtain (Appendix C) 


P.y(R) 

pv(R) 


Px(0) + 
Py( 0) - 


27t sgn cos 9 S 
1 + | cos 9 S | 
27t cos 9 S 


JY, 


1 + I cos 9 s 


JX. 


(47) 


By placing the string just above or below the equator (depending on the polarity), 
cos 9 3 = =f 0 for p = ±1, we get px = —‘ipxJY in accordance with Eq. (52) below. 
However, py = 0 is obviously wrong. 
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4-3. Conserved momenta 

To compute conserved momenta Px and Py conjugate to rigid displacements X and Y 
of the vortex with the aid of Eq. (19) for an undeformed vortex, we need the gyrotropic 
coefficients Fxy = —Fyx only (all other modes are switched off). As R = (X,Y) 
represents rigid displacements of a magnetic texture, we replace dx = —d x and dy = —d y 
to obtain the gyrotropic coefficient 

F XY = J J m • (9., : m x d y m) d 2 r = Ai tQJ (48) 

and the Poisson bracket 


{Px, Py} = -F X y = -AnQJ. (49) 

Here Q is the skyrmion charge with density 

q = m ■ (^m x 9 y m) (50) 

47T 

is its density. For a vortex with polarity p , 


Q= p/2 = ±1/2. 


(51) 


As with the domain wall, the gyrotropic coefficient for the two zero modes of a vortex 
is a topological invariant insensitive to the detailed structure of its core [25, 26]. We thus 
obtain 

P a (R) = P a (0) - 4t rQJe aP X 0 (52) 

for an undeformed vortex. Here a and (3 are Cartesian indices and e a p is the Levi-Civita 
symbol in 2 dimensions, e xy = —e yx = 1. 

Let us also compute the integral of motion associated with rotational symmetry. 
In the presence of spin-orbit coupling and dipolar interactions, rotations must involve 
magnetization components m x and m y as well as spatial coordinates x and y. To 
evaluate the conserved angular momentum of an undeformed vortex with a core cen¬ 
tered at R = {X, y), we need the gyrotropic coefficients and F$y, or F$ a for 
brevity. An infinitesimal rotation about the origin is equivalent to a global translation 
by SX a = —e a pXp A4> followed by a rotation about the vortex core. The latter operation 
does not affect the vortex because of its axial symmetry (Fig. 4). Thus a rotation reduces 
to a pure translation, giving a relation between the derivatives 


dm 9m 

9 $ =6af) a dXp 


(53) 


The gyrotropic coefficients F$x and F^y can then be expressed in terms of Fxy- 


F<s> a = ~X a F X y = -4 TrQJX a . 


(54) 


We thus obtain the angular momentum of an undeformed vortex centered at R: 


P$(R) = P*(0) 


F$ a dX a = P$(0) + 2ttQ4TR 2 ■ 


(55) 
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f.f. Conjecture of Papanicolaou and Tomaras 

Papanicolaou and Tomaras [ 6 ] conjectured that linear momentum of a two-dimensional 
ferromagnet is given by the following expression: 

P a = -4nje a p j xpqd 2 r, (56) 

Although they did not derive this result, they offered plausible arguments in favor of this 
conjecture. For example, they showed that P a is the generator of translations for the 
held m(r). Additionally, Schiitte and Garst [12] pointed out that translating a soliton by 
5X a without deformation increments linear momentum (56) by A P a = —4:TrQffe a g SXp, 
in agreement with our Eq. (52). 

Did Papanicolaou and Tomaras guess correctly the conserved momentum of a ferro¬ 
magnet in two dimensions? The answer appears to be yes. The general expression (21) 
gives the following result for linear momentum in two dimensions: 

P a (t) = -Pq(O) + J J dt J m ■ ( d a m x d t m) d 2 r. (57) 

The time derivatives dP a /dt of momenta (56) and (57) differ by a boundary term 


e a i 3 j) x f 3 m • (d 7 m x d t m) dx^. (58) 

If magnetization at the boundary is confined to a single plane (e.g., constrained to he in 
the plane of the him or to be normal to it) or if it is completely static then the boundary 
term vanishes, and the dehnition of Papanicolaou and Tomaras is equivalent to ours. 
Their conjecture is thus confirmed for reasonable boundary conditions. 

5. Discussion 

We have shown that the canonical recipe for computing conserved momenta through 
the application of Noether’s theorem quite generally fails for a ferromagnet. The problem 
is brought into focus by a close analogy between the dynamics of a ferromagnetic soliton 
and of a charged particle in a magnetic held. In the latter case, it is well known that 
canonical momenta are gauge dependent and generally not conserved. The presence of 
a background gauge held makes it necessary to follow a physical symmetry with a gauge 
transformation. Conserved momenta are generators of these combined transformations. 
They are different from canonical momenta. The group of magnetic translations [14] is 
one of the oldest examples of a gauged symmetry. It is frequently used in the context of 
the quantum Hall effect [27]. Gauged angular momentum goes even further back in time 
[28]. Wen’s projective symmetry group [29] is a relatively recent application of gauged 
symmetries. 

We have exploited this analogy to properly dehne conserved momenta of ferromag¬ 
netic solitons. Eq. (19) expresses them as a function of collective coordinates and Eq. (21) 
as a functional of the magnetization held. Computing conserved momenta P a of a soli¬ 
ton in equilibrium (no deformations) is straightforward: it requires the knowledge of 
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Figure 5: Trajectories of an interacting pair of solitons, a vortex (blue solid line) and an antivortex (red 
dashed line). Left: equal skyrmion numbers, Q\ = Q 2 - Right: opposite skyrmion numbers, Q± = —Q 2 - 


gyrotropic coefficients F a b [16, 17] of the soliton’s zero modes, which are easy to calcu¬ 
late. 

Two well-known models of ferromagnetic solitons—a domain wall in a ferromagnetic 
chain and a vortex in a thin film -have been used to illustrate these general considera¬ 
tions. We have derived the conserved momenta and have shown that canonical momenta 
obtained through the naive application of Noether’s theorem generally fail to give the 
correct answer. 

Specializing to two dimensions, our results for conserved linear and angular momenta 
confirm an earlier conjecture by Papanicolaou and Tomaras [6]. Thus our method resolves 
paradoxes accumulated over decades and provides a straightforward and very general 
method for computing conserved momenta of ferromagnetic solitons. 

It is worth noting that early works by Slonczewski [1] and by Thiele [2] contained 
the correct treatment of linear momentum of a domain wall. These researchers relied 
on Newton’s second law for solitons and calculated linear momentum by integrating the 
impulse of an external force along the lines of Eqs. (6) and (7). They even anticipated 
“the treacherous nature” [2] of canonical momentum (44)! 

The existence of conserved momenta can be useful in situations where solitons interact 
with one another. For example, if two vortices in a thin film (typically a vortex and an 
antivortex) attract each other, their net conserved linear momentum is 

P = -Anjz x (QiR.! + Q 2 R 2 ). (59) 

If the vortices have equal skyrmion numbers, Q\ = Q 2 , then their average position 
remains in place, (Ri + R 2 )/2 = const, which means that they orbit a common center. 
If the skyrmion numbers are equal and opposite, Q\ = —Q 2l then their relative position 
remains fixed, Ri — R 2 = const, as the two vortices move in the direction orthogonal 
to the line connecting them. The two situations are illustrated in Fig. 5. The gradual 
convergence of the two solitons in the figure reflects the influence of weak viscous friction. 

Similarly, two consecutive domain walls in a ferromagnetic wire (Sec. 3) have the 
total linear momentum Pz = 2J'(4?i — d> 2 ). In the absence of external forces, the 
two interacting domain walls will precess at the same frequency but keep their relative 
orientation fixed, $1 — $2 = const. 
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Appendix A. Magnetic field F^ from gauge potential Aj 

The magnetic field F i? felt by a soliton (20) can be derived in a number of ways. 
Tretiakov et al. [16] obtain it by translating the Landau-Lifshitz equation of motion for 
the magnetization field m(r,i) into the language of collective coordinates q (t). Here we 
obtain it directly from the definition of the gauge potential (8) and magnetic field (9): 


F — 

-*■ in — 


dAj dAi f da a dm a 

dq,, dqj J 

da a dmp dm c 
dmp dqt dqj 

fa-, „ dm P dm < 

J J a/37 / dq x dqj 


dqi dqj 




dV = 


dV = J 


dV 

/ da a _ dap \ dmp dm a 
\ dmp dm a J dqi dqj 


m 


/ dm dm 
\ dqi dqj 


dV. 


(A.l) 


In this derivation we relied on Eq. (1), which states that the spin vector potential a(m) 
describes a magnetic monopole, V m x a = — J m, or 


da a dap 
dmp dm a 


\J ^a.p'i'^a^f • 


(A.2) 


Appendix B. Poisson brackets for collective coordinates and conserved mo¬ 
menta of a soliton 


We first establish that the tensor of Poisson brackets for collective coordinates {qi, qj} = 
Qij is the inverse of the gyrotropic tensor (20): FQ = 1, i.e., FijQjk = Sik- (Sum¬ 
mation over doubly repeated indices is implied.) 


Fij Q jk 


dVdV' 


f dcj) d cos 9 


\dqi 


dV 

K UQi 0(j) 


dqj 

dcos 9 


dqi 


d cos 9 d(f> \ | 
dqi dqj J ' 
Sqk \ _ dqk 
6 cos 9 ) dqi 


dqj 6q k 
S cos 9' 5(j)' 

= dik- 


dqj dq k \ 
dcj) 1 dcos9’) 


(B.l) 


Here we used the differentiation chain rule and functional derivatives 


dcos 9 
d cos 9’ 


dcj) 

dfjy 


<5(r - r'), 


dcos 9 

d# 


dcj) 

d cos 9' 


(B.2) 


We thus obtain the Poisson bracket for two collective coordinates in Eq. (27), 

{q i ,q j } = (F- 1 ) ij . (B.3) 
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The remaining two results in Eq. (27) can be obtained by using the identity 


tl9} = ^,W, {q " q,) (B4) 

valid for arbitrary functions /(q) and g( q). The mixed Poisson bracket in Eq. (27) is 

{Pa,Qi} = = - F a j6 ik (F ~ 1 ) jk = - Sai , (B.5) 

OQj UQk 


where the identity dP a /dqj = — F a j follows from the definition of conserved momentum 
(19). 

The Poisson bracket for two conserved momenta in Eq. (27) is obtained along the 
same lines: 


{Pa,Pb} 


dPadPb 

dq t dqj 


{ Qi ; Qj } — PaiPbj(P )ij 


—Fab- 


(B.6) 


Appendix C. Canonical momentum of a vortex 

We evaluate the canonical momentum p of a vortex, Eq. (44). The integrand decays 
slowly away from the core: cos 0± 1 — > ±1 and d a (f> ~ —e-apxp/r 2 . To ensure convergence, 
we restrict the integration to a finite disk of radius Rd greatly exceeding both the core 
size and displacement R: 

p Q (R) = — J I (cost? ± l)d a cj)d 2 r. (C.l) 

J r<R d 

Next we subtract the momentum of a vortex with the core at the center of the disk, 
whose configuration is m 0 (r): 


p Q (R)-p a (0) = -J / (cos0 o ± l)d a <f> 0 d 2 r 

J\r-R\<R d 

+J f (cos 0 O ± l)9 a ^o d 2 r. (C.2) 

Jr<R d 

The contributions from the area where the disks r < Rd and |r — R| < Rd overlap cancel 
out and we are left with an integral over two crescents of thickness R bounded by the 
edges of the disks. For a disk radius Rd greatly exceeding the core size, we may set 
cos 0o 0 in the crescents. 

The integration over r is conveniently done in polar coordinates, radius r and direction 
n = r/r = (cos sin ip). For Rd R , d a (f>o = e a pnp/r is approximately constant across 
the thickness of the crescent, so that for a given n, the integration over r yields a factor 
of n • R = n 7 X 7 . Integration over directions yields Eq. (46): 

Pa{ R) - Pa(0) = ± J / e a pnpn 1 X 1 d'tjj = ±irJe a 0 Xp, m s = ±z. (C.3) 

Lastly, we derive the canonical linear momentum when the Dirac string is attached at 
a generic location along the (j) s = 0 meridian, m s = (sin0 s , 0, cos0 S ), with the gauge (3). 
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We use the same approach as above, introducing a long-distance cutoff Rd and evaluating 
the momentum difference p(R) — p(0). Doing so again reduces the integration area to 
two narrow crescents of radius Rd and width R. Since the integration area is far away 
from the core, magnetization lies in the xy plane and 


m s • (mo x <9 Q m 0 ) = m s • zd a (po = - cos O s e a pnp/R d . 
We thus obtain Eq. (47): 


P.y(R)- px(O) = 


Py(R)-Py(0) = 


_ f 2n cos 9 S sin ip(X cos ip + Y sin ip) 

+J / - : —:—;— 7 —:—7 - dip 

Jo 1 + sin 9 S sin ip 

2ir sgn cos 9 S 

1 + | cos 6 s | J ’ 

_ /' 27r cos 9 S cos ip(X cos ip + Y sin ip) 

-J I -——. „—t —7 - dip 


J 0 

27t cos 9 s 
1 + I cos 9 


1 + sin 9 S sin ip 


rJX. 


(C.4) 


(C.5) 
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